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The branching fraction of the tp" decays into charmless final states is estimated in the S- and 
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states. The experimental search for these decays is also discussed. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.40.Gx 



I. INTRODUCTION 

As more and more data are collected at ?/'(3770) (short- 
ened as 27.7 pb-i by BES-II, 55 pb^^ by CLEOc, 
and a few fb^^ is in progress by CLEOc, the study of the 
non-DD decays of V'" gets renaissance recently. 

Due to its mass coincides with the DD mass threshold, 
-0" is believed to decay predominantly to DD, and the 
non-DD decays including transitions between charmo- 
nium states and decays into light hadrons are expected 
to be small. The first exclusive non-DD decay mode of 
ip" was searched for using the Mark-Ill data sample 1], 
where evidence of ip" — > J/i/^tt+tt" was observed; after- 
wards theoretical calculation of its decay width was pub- 
lished 0. This decay mode was searched for recently 
by BES and a signal of 3.1a was observed Q, which 
is in marginally agreement with the upper limit set by 
CLEOc The other searches of the non-DD decays 
of ^" were all reported in either doctoral theses [l], Q 
based on the Mark-Ill data or conference talk based 
on the BES-II data, and no statistically significant re- 
sults were given. This may indicate that the data sam- 
ples utilized are still not large enough to search for these 
channels of small branching fractions. Now with much 
larger -0" samples, measurements of the cross section 
a{e~^e~ ip" DD) = a{DD) have recently been 
reported by BES-II 0] and CLEOc collaborations us- 
ing either single-tag or double-tag method. Comparing 
with the measured cr(e"*'e^ — > V" ~^ anything), the dif- 
ference can be up to 1.4 nb (about 18% of the total cross 
section of the ip" production in e+e" collision 0), which 
indicates the existence of substantial non-DD decays for 
-0". However, there is limited theoretical work on this 
aspect. In Ref. it is estimated that at most 600 keV 
(- 2.5%) of the V'" total width of (23.6 ±2.7) MeV is due 
to the radiative transition, and perhaps as much as an- 
other 100 keV (~ 0.4%) is due to the hadronic transition 
to J/tpTTTT. All these together are far from accounting for 
a deficit of 18% of the total ip" width. 

In this paper, we concentrate on the charmless decays 
of tp". By charmless decay, we exclude those decay modes 
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with either open or hidden charm. We shall estimate 
the partial width of the charmless decay of ip" in the 
2S'-1D charmonia mixing scenario, under the assumption 
that the perturbative QCD (pQCD) "12% rule" holds be- 
tween pure IS* and 25 charmonium states. The mixing 
between 2S and ID states of charmonium was originally 
proposed to explain the measured large Fee of ip" ■ It is 
suggested that the mass eigenstates ■0(3686) (short- 
ened as ip') and ip" are the mixtures of the 25"- and ID- 
wave of charmonia, namely '0(2^ 5*1 ) and ■ip{l^Di) states. 
This was used by Rosner [l3| in explaining the "p7r puz- 
zle" in tp' and J/0 decays. He suggested that the mixing 
of ■0(2'^S'i) and 0(l'^il'i) states is in such a way which 
leads to almost complete cancellation of the decay am- 
plitude of -0' — > pTT, and the missing pir decay mode of -0' 
shows up instead as enhanced decay mode of tp" . This 
idea was then applied to solve the enhanced decay of 
tp' KgK^ relative to J/tp KgK^, and a prediction 
of the branching fraction of 0" KgR^ was given l l}. 
In principle, if this scenario is correct, it can be gener- 
alized and relate the partial widths of each individual 
mode in J/ip, ip' and tp" decays. Therefore by virtue of 
the branching fractions of J/0 and tp' decays, we are able 
to estimate the corresponding decay branching fraction 
of tp". This will be tested by the upcoming large tp" data 
sample from CLEOc. 

In the following parts of the paper, we begin our study 
with a general review of the 12% rule between tp' and 
J/tp decays, followed by an introduction of the 25'- and 
ID- wave charmonia mixing scheme. Then we conduct 
the calculation on the possible charmless decay width of 
tp" with the available information from J/tp and tp' de- 
cays. We show that the scenario and currently available 
information on J/tp and tp' decays can accommodate a 
partial width of 3.0 MeV charmless decays of tp". Finally 
we turn to the experimental searches for these decays. 
Both inclusive and exclusive methods are examined. 



II. 12% RULE AND 2S-ID MIXING 

From the pQCD, it is expected that both J/tp and tp' 
decaying into light hadrons are dominated by the annihi- 
lation of cc into three gluons, with widths proportional to 
the square of the wave function at the origin |\E'(0)p [T^ . 
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This yields the pQCD 12% rule, that is 



Qh = 



B 



B 



'J/ ip — >-e+e 



12.7%. 



(1) 



The violation of the above rule was first observed in 
pTT and K*^K^ +c.c. modes by Mark-II JJ|j. Since then 
BES-I, BES-II and CLEOc have measured many two- 
body decay modes of V'' El IH III IH El 111 HIM 111- 

Among them, some obey the rule , like baryon-antibaryon 
{BB) modes, while others are either suppressed in ip' 
decays, like vector-pseudoscalar (VP) and vector-tensor 
(VT) modes, or enhanced, like KgKl. There have been 
many theoretical efforts trying to solve the puzzle 10, 231. 
Among them, the 25- ID charmonia mixing scenario (10l | 
predicts with little uncertainty B{ip" — > pvr) which agrees 
with experimental data [ij, |2J| . 

In S- and D-wave charmonia mixing scheme, the mass 
eigenstates ip' and ip" are the mixtures of the S- and 
D-wsLve charmonia, namely 

IV'') = |235i)cos6i- |l3i:»i)sin6' , 
'-I- |l3i:)i)cos6i , 



smf 



where 9 is the mixing angle between pure ipi'^^Si) and 
-0(1^-01) states and is fitted from the leptonic widths of 
V-" and V' to be either (-27 ±2)° or (12 ±2)° T^. The 
latter value of 6 is consistent with the coupled channel 
estimates H 12^ as well as the ratio between ip' and ip" 
partial widths to J/V'tt+tt" 0- Hereafter, the calcu- 
lations and discussions in this paper are solely for the 
mixing angle 6* = 12° |2^ . 

As in the discussion of Ref. E3j since both hadronic 
and leptonic decay rates are proportional to the square 
of the wave function at the origin, it is expected that if 
ip' is a pure ip{2'^Si) state, then for any hadronic final 
states /, 



ri^' ^ /) = r(j/v. ^ /) 



r(0' ^ e+e-) 
T{J/^P e+e-) 



(2) 



The electronic partial width of J/iJj is expressed in po- 
tential model by [23 

r(J/^_e+e-) = ^|i?is(0)r 

J/i> 

with a the QED fine structure constant, Cc = 2/3, Mj/^ 
the J/4> mass and i?is(0) the radial l^Si wave function 
at the origin. 

Since tp' is not a pure ip(2^Si) state, its electronic par- 
tial width is expressed as [W^ 

Wei 



r(^'^e+e-) = 



Ml, 



cos6'i?2s(0) 



2^/2r, 



■ sin0R'ljj{O) 



with M^' the mass, rric the c-quark mass, i?2s(0) 
the radial 2^ Si wave function at the origin and R'(jj{0) 



the second derivative of the radial l^Di wave function 
at the origin. In the calculations in this paper, we 
take i?2s(0) = 0.734 GeV^/^ and 5i?'/^(0)/(2V^m2) = 
0.095 GeV3/2 from Refs. Ellli- 

If Eq. holds for a pure 2S state, j^' ^ f, ^' ~^ f 
and J/ip — > / partial widths are to be 



C 

|sin6ii?2s(0) +7ycosf'i^ 



nr - 
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r(V'- 


-/) 


r( j/^ - 





C 



^\coseR2s{O)-r]sm0f , 



Ml, 



\Rism' 



(3) 



where C/ is a common factor for the final state /, M^,, 
the ip" mass, and ry = |ry|e"^ is a complex parameter 
with (f> being the relative phase between {f\\^Di) and 
(/|2^5i). 

From Eq. it is obvious that with r(J/i/' — + /) 
and T{il>' — > /) known, two of the three parameters, C/ 
and complex 77, can be fixed, thus can be used to predict 
r('0" — > /) with only one unknown parameter, say, the 
phase of 77. Thus the S- and D-wave mixing scenario 
provides a mathematical scheme to calculate the partial 
width of 4'" decay to any exclusive final state, with its 
measured partial widths in J/tp and ip' decays. 

However, the current information concerning the ip' 
decay is extremely limited, which prevents us from esti- 
mating Qh values for most exclusive decay modes. Ta- 
ble lists some hadronic final states which are measured 
both in J/ip and tp' decays, together with the calculated 
Qh defined in Eq. (Q. Summing up all the channels in 
Tabled makes less than 2% of the ip' decay through ggg 
annihilation. 

From Tabled we notice that compared with the 12% 
rule, the -0' decays to 

1. the pseudoscalar-pseudoscalar (PP) mode KgKl is 
enhanced; 

2. the VP and VT modes are suppressed; 

3. most of the BB modes are consistent with it. 

The summed branching fractions and Qh values for these 
three categories of decay modes are evaluated and also 
listed in Table In estimating the charmless decays of 
ip", we shall discuss these three different cases separately. 

Since the experimental information on the exclusive 
decays of V' is rather limited, we turn to inclusive 
branching fractions of J/ip and tp' hadronic decays as 
an alternative. The estimation is based on the assump- 
tion that the decays of J/tp and ip' in the lowest or- 
der of QCD are classified into hadronic decays {ggg), 
electromagnetic decays (7*), radiative decays into light 
hadrons {jgg), and transition to lower mass charmo- 
nium states {ccX) |^ Is^l- Thus, using the relation 
B{ggg) + Bi-fgg) + B{'j*) + B{ccX) = 1, one can derive 
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TABLE I: Branching fractions and Qh values for some J /ij) and '4>{2S) decay channels. 



Modes 


Channels 


Sj/^(10-^) 


S^(2S)(10~*) 


Qh (%) 


Rcf. 


0^0" 


K+K- 


0.147 ± 0.023 
0.237 ±0.031 
0.182 ± 0.014 


0.8 ± 0.5 
1.0 ±0.7 
0.52 ±0.07 


54 ± 35 
42 ±30 
28.8 ±3.7 


[29] 

M 

[16, iU 


sum 


PP 


0.57 ±0.07 


2.32 ± 1.27 


41.1 ± 22.8 




1~0~ 


K^'K* {m2f + c.c. 


12.7 ± 0.9 
5.0 ± 0.4 

4.2 ± 0.4 
0.42 ±0.06 


0.29 ± 0.07 
fl 1 5 ± n 08 
1.10 ±0.20 
0.20 ±0.06 


0.23 ± 0.6 
0.3 ± 0.2 
2.6 ±0.5 
4.8 ± 1.6 


[301 


1-2+ 


c^/2(1270) 
pa2 

K* (892)° K*2{U30)° + c.c. 
<?i/2(1525) 


4.3 ±0.6 
10.9 ±2.2 
6.7 ±2.6 
1.23 ±0.21 


2.05 ±0.56 
2.55 ±0.87 
1.86 ±0.54 
0.44 ±0.16 


4.8 ± 1.5 
2.3 ± 1.0 
2.8 ± 1.3 
3.6 ± 1.4 


L15J 


sum 


VP& VT 


45.45 ± 7.37 


8.64 ±2.54 


1.90 ±0.64 




BB 


PP 
AA 

E(1385)=^E(1385)+ 
A++A~- 


2.12 ± 0.10 
1.30 ± 0.12 
1.27 ±0.17 
1.03 ±0.13 

1.8 ±0.4 
1.10 ±0.29 


2.07 ± 0.31 
1.81 ± 0.34 
1.2 ±0.6 
1.1 ±0.4 
1.88±0.62t 
1.28 ±0.35 


9.8 ± 1.5 
13.9 ± 2.9 
9.4 ±4.9 
10.7 ±4.1 
10.4 ±4.2 
11.6 ±4.4 


[22J 


sum 


BB 


8.62 ± 1.21 


9.34 ±2.62 


10.8 ± 3.4 






Note: t simple 


normalization by 


' = (1/2)HH. 







^{999) + ^{199) by subtracting B["f*) and B{ccX) from 
unity. 

The calculated values of .6(7*) and B{ccX)^ together 
with the values used to calculate them are summarized 
in Table im As regards to -0', two final states 777(26') and 
hc{l^Pi)+X with faint branching fractions are neglected 
in our calculation. By deducting the contributions B{'^*) 
and B{ccX), we find that 6(J/V' ^ 999) ± B{J /ii 
-igg) = (73.5 ± 0.6)% and B{^' ^ ggg) + S(V/ ^ 755) = 
(19.1 ± 2.5)%, then the ratio of them is 



Q9 



Bjip' ggg + igg) 
B{J/ip ggg + -fgg) 



(26.0 ±3.5)%. (4) 



The above estimation is consistent with the previous 
ones The relation between the decay rates of 

ggg and 7017 is readily calculated in pQCD to the first 
order as ISal 
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T{J/^p jgg) ^ 

T{J/^ ggg) 5 ^^(toc) 



fl-2.9^ 



Using as{mc) = 0.28, one can estimate the ratio to be 
0.062. A similar relation can be deduced for the ijj' de- 
cays. So we obtain B{J/^ — > ggg) ~ (69.2 ± 0.6)% and 
B{i^' ggg) (18.0 ±2.4)%, while the "26.0% ratio" in 
Eq. stands well for both ggg and jgg. Although Qg 
is considerably enhanced relative to Qh in Eq. , it co- 
incides with the ratio for the RgK^ decay mode between 
ip' and J/ipj which is 



according to the recent results from BES 

MM The 

relation in Eq. Q was discussed in the literature as the 
hadronic excess in ip' decay [s^, 13 ■ It implicates that 
while some modes are suppressed in ^p' decays, the dom- 
inant part of Ip' through ggg decays is enhanced relative 
to the 12% rule prediction in the hght of J/ip decays. 



TABLE II: Experimental data on the branching fractions for 
J/tp and tp' decays through virtual photon and to lower mass 
charmonium states used in this analysis. Most of the data are 
taken from PDG (2^, except for B{J/ip, ip' ^ "f* ^ hadrons), 
which are calculated by the product R ■ B{J/tp,tp' — » ^'^ fi~), 
with R = 2.28 ± 0.04 [sl. In estimating the errors of the 
sums, the correlations between the channels are considered. 



Channel 


s(J/V) 


S(V') 


7* hadrons 


(13.4±0.33)% 


(1.66±0.18)% 


e+e" 


(5.93±0.10)% 


(7.55±0.31)xl0-^ 




(5.88±0.10)% 


(7.3 ±0.8) X 10"^ 


T+T- 




(2.8 ±0.7) X 10"=' 


j* ^ X 


(25.22±0.43)% 


(3.43±0.27)% 


yvc 


(1.3±0.4)% 


(2.8±0.6)xl0"^ 


-k'^-k~ J/tp 




(31.7±1.1)% 






(18.8±1.2)% 


■qj/ip 




(3.16±0.22)% 


-k'^J/iP 




(9.6±2.1)xl0"'* 


7XcO 




(8.6±0.7)% 


7Xci 




(8.4±0.8)% 


7Xc2 




(6.4±0.6)% 


ccX 


(1.3±0.4)% 


(77.4±2.5)% 



Qj,OKO =(28.8 ±3.7)% 



(5) 
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III. THE CHARMLESS DECAYS OF 7/." 

We define the enhancement or suppression factor as 



Qif) 



(6) 



In the 2S-1D mixing scheme, for any final state /, its 
partial width in ip" decay can be related to its partial 
width in J/ip and ip' decays by Eq.®, with an unknown 
parameter which is the phase of 77. This unknown phase 
constrains the predicted ^ /) in a finite range. We 
calculate 



Rr ^ r{r - /)/r( J/v - /) 



(7) 



as a function of Q{f) and plot it in Fig. ^ In the figure, 
the solid contour corresponds to the solution with (p — 0; 
the dashed one corresponds to the solution with (p = 
180°; and the hatched area corresponds to the solutions 
with (j) taking other non-zero values. 

To make it clear, we discuss the final states in three 
situations: Q{f) < 1, Q{f) > 1, and Q{f) = 1. 



A. Final states with Q{f) < 1 

If Qif) < the decay ^p' ^ f is suppressed relative 
to J/V' /. The extreme situation is Q{f) 0, cor- 
responding to the absence of the mode / in ip' decays. 
This is the case which was assumed for the pn mode in 
the original work to solve the pn puzzle by the S- and D- 
wave mixing If Q(/) = Oi the solution of the second 
equality of Eq. ^ simply yields ry = _R2s(O)cos0/ sm9 
which cannot have a non-zero phase, and Rr — 9.2. 

Generally, the suppression factor could be different 
from zero. Even the pn and other strongly suppressed 
VP modes are found in -0' decays recently by BESII 
and CLEOc ^ with QiVP) ~ 0(10-^). In this case, 
there are two real and positive solutions of rj as shown in 
Fig. n corresponding to the maximum and minimum of 
their possible partial widths in ip" decays. The solutions 
with 77 having a non-zero phase yield the values of Rr 
between the minimum and maximum limits. 

For VT final states, which are measured to have 
QiVT) « 1/3, with Eq. ®, we get 2.0 < i?r < 21.6, 
as shown in Fig. ^ where the upper and lower limits cor- 
respond to two real and positive solutions of rj, the range 
is due to the values of 77 with non-zero phases. 



B. Final states with Q{f) > 1 

If Qif ) > Ij the decay ■0' — > / is enhanced relative to 
J/ip f. The extreme situation is Qif) — > 00, corre- 
sponding to the complete absence of the final state / in 




FIG. 1: -Rr versus Q{f). The solid contour corresponds to 
<?!> = 0; the dashed contour corresponds to <p = 180°; and the 
hatched area corresponds to (p having other non-zero values. 
The inset displays the variation of Rr in the vicinity of Q(/) = 
1, where the dot-dashed line denotes i?r = 1. 



J/^ decays. For 

Qif)> 



cos6'i?2s(0) 



cos0i?2s(O)-^^sin^?i?'/^(O) 



= 1.06, 



there are two real solutions of "q, one is positive and the 
other negative. From the first equality of Eq. it is 
seen that the positive solution leads to the larger Rr- i.e. 
larger r(-0" f) (the solid contour in Fig. while the 
negative solution leads to the smaller Tiip" —* /) (the 
dashed contour in Fig. P). 

For the known enhanced mode in ip' decays like KgR^, 
Q{K°Kl) = 2.26, we find 1.4 < i?r < 52.5, which cor- 
responds to the -0" decay partial width from 0.024 to 
0.87 keV. 

It should be noted that for finite r(V/ /), Qif) ^ 
00 means the diminishing of TiJ/tp — > /), which gives 
rise to Rr 00 according to the definition of Rr in 
Eq.©. Under such circumstance, it is more intuitive to 
calculate 

R^ ^ Fir ^ /)/r(7/^' ^ /). 

Its variation as a function of Qif) is shown in Fig.|21 As 
Qif) — > 00, the solid and dashed contours converge into 
the same point (i?p — > 21). In such case S-wave state 
does not decay to /, but D-wave does. Its partial width 
in ip' and ^" decays comes solely from the contribution 
of 77, or the D-wave matrix element. 



C. Final states with Q(/) = 1 

These final states observe the 12% rule in J/tp and 
Ip' decays. In this case, apparently one solution is 
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FIG. 2: 7?p versus Q{f). The solid contour corresponds to 
(j) = 0; dashed contour corresponds to </> = 180°; and the 
hatched area corresponds to (j) having other non-zero values. 



identical to the pure electromagnetic decays, with r/ ~ 
5R'ljj{0)/{2V2ml) and i?r = 0.048. As in the leptonic 
decays, their partial widths in ip" decays are small rela- 
tive to the partial widths in J/ip decays. However, there 
are also other solutions with overwhelming contribution 
from 77 which lead to very large partial widths in ■0" de- 
cays. This can be seen from Fig. ^ at the point with 
Q{f) = 1. We get the maximum of Rr of 34.0 which 
corresponds to another real and positive value of rj. If rj 
has a non-zero phase, then 0.048 < Rr < 34.0. 



D. Numerical results 

From Fig. ^ we see that except in the range 0.52 < 
Q{f) < 2.06 and a small range of the phase, Rr is always 
greater than 1. This range excludes virtually all known 
decay modes except BB which has Q{BB) « 1. Even 
inside this range, there are other solutions by which Rr 
(or i?p) is at 0(10). It means in this scenario, contrary 
to the naive guess, the charmless decay width of -0" is 
greater than that of J/-0 or . More surprising is that 
Rr could be as large as a few tens for the final states with 
Q{f) > 1. In general, the final states which are enhanced 
in ip' decays possibly have a large combined partial width 
in ip" decays, especially if the phase of rj is zero or very 
small. 

There are reasons to assume that the phase 4> which is 
between the matrix elements (/|2S') and {f\lD) should 
be small 11]. For the decay mode like pn, since there 
is almost complete cancellation between cos6'i?2s(0) and 
■qsu\6 so that {p-K\tp') = cos^i?2s(0) — -qsinO « 0, the 
phase of 77 must be small. If this is to be extrapolated to 
all final states, the physics solution will follow the solid 
contour of Fig. ^ Another argument comes from the 



universal phase between the strong and electromagnetic 
amplitudes of the charmonium decays. It has been known 
that in the two-body decays of J/V'j the phase between 
the strong and electromagnetic amplitudes is universally 
around 90° |3lll33| . Recently, it has been found that this 
phase is also consistent with the experimental data of ^' 
and V" decays 123, [3^. Since there is no extra phase 
between 25* and ID matrix elements due to electromag- 
netic interaction, as in the calculations of the leptonic 
decay rates of ij}' and V'": the universal phase between 
the strong and electromagnetic interactions implies there 
is no extra phase between the two matrix elements due 
to the strong interaction too, i.e. (/) « 0. This conclusion 
means, for the modes which are enhanced in -0' decays, 
their partial widths in ip" decay must be greater than 
those in J /tp or -0' decays by more than an order of mag- 
nitude. 

In Sect.|nl we estimate that B{J/4> — > ggg) — (69.2 ± 
0.6)% while 6(0' -> ggg) ~ (18.0 ± 2.4)%. Among 
the final states, we know that VP and VT modes have 
Q(/) < 1. For them, 

ES(J/i/' ^ VP, VT) « 4.6% , 
^6(V'' ^ VP, VT) « 8.6 X lO-'' . 

Furthermore, there are final states with Q{f) ~ 1 (such 
as BB), for them, 

Y.B{J/iP^ BB) w 0.9% , 
X;S(0' ^ BB) « 9.3 X 10-4 . 

After subtracting the final states which are known to have 
Q{f) < 1 and Q{f) w 1, the remaining 63.8% of Jf-ip 
decay with a total width of r(J/-0 r.f.s.) « 58.1 keV, 
and 17.8% of 0' decays with a total width of r(0' 
r.f.s.) « 50.1 keV which are gluonic either has Q{f) > 1 
or Q{.f) unknown. Here r.f.s. stands for the remaining 
final states. On the average, these final states have 

Qir.f.s.) « 2.19 . 

This is roughly comparable to QiK^Kl) = 2.26. 

If there is no extra phase between 25 and ID matrix 
elements as argued above, then Rr takes the maximum 
possible value with Rr ~ 51.6 for Q{r.f.s.) « 2.19. With 
this value, we find r(0" ^ r.f.s.) = Rr x T{J/i; 
r.f.s.) « 3.0 MeV for the partial width of these remaining 
final states in ip" decays, which is 13% of the total 0" 
width. 

The above calculation takes the averaged Q{f) for the 
final states with Q{f) > 1 and Q{f) unknown, so it 
merely serves as a rough estimation. The exact value 
of the partial width should be the sum of the individual 
final states which in general have various Q{f) values. 
At present a major impediment to do the accurate evalu- 
ation is the lack of experimental information. Neverthe- 
less, if we take 13% as charmless decays (the calculations 
done channel by channel for VP, VT and BB modes in 
Table give a summed maximum possible width of 93 
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keV in ip" decays, or 0.4% of the V'" total width) together 
with the charmonium transition contributions of 3% 
(2.5% for radiative transition and 0.4% for mrJ/tp), we 
obtain a maximum of ijj" non-DD decay branching frac- 
tion of 16% in the 2S-1D mixing scenario to be compared 
with 18% as summarized in Ref. Hi. 



IV. EXPERIMENTAL MEASUREMENTS 

With the expected data of a few fb~^ at V'" from the 
running CLEOc and more from the future BES-III, it is 
of great importance to search for the charmless decays 
of V'" experimentally. These measurements test the 25- 
ID mixing scenario on the explanation of the pir puzzle 
in J/^p and ip' decays, and provide information on the 
relative phase between the 25 and ID matrix elements, 
which is related to the charmonium decay dynamics. 

Besides the theoretical interests, such measurements 
are important for the experiment itself. So far, in the 
fitting of the ip" resonance parameters with the scanned 
cross sections, it has been assumed that ^p" decays com- 
pletely into DD [Hill m 113 . If the non-DD branching 
fraction is substantial, it must be considered in the fitting 
in order to get self-consistent results. 
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FIG. 3: The cross section a{e^e^ hadrons) in the vicinity 
of the tp" resonance calculated with parameters provided by 
PDG. The total cross section (t*°' is divided into four parts: 
the non-resonance aMR, the radiative tails of J/tp {oj/^) and 
V'' i'^ip'), and the ip" resonance (cr^"). 



A. Inclusive method 

A substantial non-DD decays of -0", including charm- 
less decays and charmonium transitions, originally 
caught attention from the comparison between the cross 
sections of the inclusive hadrons and DD at the ■0" peak. 
However, this is not unambiguous due to the poor statis- 
tics of the data samples and the complexity of the anal- 
ysis. In addition, results from different experiments are 
consistent with each other only marginally. 

Fig-Elshows a{e^e~ hadrons) in the vicinity of the 
■0" resonance calculated with parameters by PDG [2^. 
The total cross section (T*°* can be expressed as 



(8) 



which contains four parts: the non-resonance cross sec- 
tion a MR, the radiative tails of J/ip (cj/^) and ip' (cr^')i 
and the ip" resonance (cr^"). The non-resonance cross 
section is usually expressed in terms of R value and the 
/i pair cross section as gmr — R-u{e^e~ — > fi~^ii~). The 
Breit-Wigner formula is adopted to depict resonances of 
J/ip, ip' and ip" , where the total decay width of ip" is 
energy dependent: 



2 

cm 



with 



r^"(£'cm) — Cy 



Ml„) 



Pdo 



tl)" 



■n3 



1 + {rpjjoY 



l + (rp^±) 



(9) 



where p is the D'^ or D^ momentum, r is the classical 
interaction radius, and Cr is defined as follows: 



Cv 



Pdo 



Pd± 



=M,,, 



Here F,;,// (M,;// ) is the ip" total decay width given by 
PDG [23. 

In previous analyses |37ll38Ll39ll40| . the observed inclu- 
sive hadronic cross sections were fit with the theoretical 
one with the contributing terms in Eq. ((SJ to obtain the 
resonance parameters, which yield the ip" cross section 
at the peak. The comparison of this cross section with 
the DD cross section measured by tagging D mesons 
yields the non-DD decay branching fraction. An incon- 
sistency exists in this procedure because in the fit to the 
total cross section, DD final states was assumed to satu- 
rate the ip" decays. Since light hadrons have much lower 
thresholds than DD, a large non-DD branching fraction 
directly affects the shape of the resonance curve, and also 
indirectly through the energy-dependent width. Taking 
into account the non-DD decays, Eq. lO should be re- 
vised by including another term, that is 



r^"(£'cm) — Cp X 



Ppa 



pI± 



non-DD 



l + irp^oY l + irp^±y 
where C'^q^.^^-jj is proportional to the part of the width 
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from non-DD decays, and 



non-DD 



With the Cnon-Z5"D term in the expression for T^n , the 
fitting of the resonance curve to extract the resonance 
parameters must be done together with the fitting of the 
DD cross section. In this procedure, the non-DD decay 
branching fraction is extracted together with the reso- 
nance parameters. 

However, this method subtracts DD cross section from 
the total inclusive one to get the non-DD cross section 
which is only a fraction of the total inclusive cross section. 
It suffers from the large uncertainties of the measure- 
ments, so it is difficult to obtain a statistically significant 
result. 



B. Exclusive method 



The calculations in section IIIII show that those final 
states which are suppressed in ip' decays relative to J/tp, 
and especially those ones enhanced in ijj' decays, will 
show up in ip" decays with maximum possible partial 
widths more than an order of magnitude greater than 
their widths in J/ip or ip' decays. So these exclusive 
charmless modes should be searched for in -0" decays. 
This provides direct test of the calculations based on the 
2S-1D mixing scheme. Here we discuss three typical ex- 
clusive modes: VP mode which is suppressed in ?/;' de- 
cays relative to J/V', PP mode which is enhanced, and 
of particular interest, the BB and (/)/o(980) modes which 
observe the 12% rule. 

To measure the exclusive VP mode in ^" decays by 
e~^e~ experiments, the contribution from non-resonance 
virtual photon amplitude and its interference with the 
resonance must be treated with care. A recent study 
on the measurement of tjj" — s- VP in e^e~ experiments 
shows l2j| that with the decay rate predicted by the S- 
and D-wave mixing, the interference between the three- 
gluon decay amplitude and the continuum one-photon 
amplitude leads to very small cross sections for some VP 
modes, e.g. pir and K*^K~ -\-c.c., due to the destructive 
interference, but much larger cross sections for other VP 
modes, e.g. K*°K -\-c.c. due to the constructive interfer- 
ence. In another word, although the branching fractions 

of p°7r" and if ° differ by only a fraction due to SU(3) 
symmetry breaking |4l| . their production cross sections 
in e~^e~ collision differ by one to two orders of magnitude. 

Among the PP modes, there is the KgR^ final state 
which decays only through strong interaction and does 



not couple to virtual photon |llll41| . There is no compli- 
cation of electromagnetic interaction and the interference 
between it and the resonance. So the observed KgKl in 
e~^e~ experiment is completely from resonance decays. 
In the 2S-1D mixing scheme, with the BES recently re- 
ported KgKl branching fractions in J ftp and tp' [iTj 
decays as inputs, it is estimated ^ (1.2 ±0.7) x 10~^ < 
B{ip" K^Kl) < (3.8 ± 1.1) X 10-5 113. If there is no 
extra phase between {K^Kl\2^Si) and {K^Kl\l^Di), 
then its branching fraction is at the upper bound. With 
17.7 pb"^ V'" data, BES has set an upper limit jl^l, which 
is still beyond the sensitivity for testing the above pre- 
diction. More precise determination of this branching 
fraction is expected from the analysis based on larger 
data samples of CLEOc and BES-III. 

However, for other PP modes, or more generally other 
final states which are enhanced in Tpi' decays, in e^e^ 
experiments there is still the complication from the non- 
resonance virtual photon amplitude and its interference 
with the resonance. 

Of particular interest are the final states with Q{f) ~ 
1. These are the BB modes and the vector-scalar mode 
Q. As discussed in subsection rill Cl for 0( f] — 



1, there are two real and positive solutions with i?r = 
0.048 and i?r = 34.0. These two solutions are three 
orders of magnitude apart, their branching fractions are 
extremely sensitive to the relative phase between the 25* 
and ID matrix elements. The BB branching fractions 
in J/-0 decays are at O{10'^), while 0/o(98O) is (3.2 ± 
0.9) X IQ-'*. If the physics solution of 77 is the larger one 
of the two real values, then the BB branching fraction 
in ip" decay would be at ©(lO""*) and 0/o(98O) would be 
4.2 X 10~5, which can be observed in the ip" data sample 
over 1 fb^^. 



V. SUMMARY 

Based on the available experimental information of 
J/ip and ip' decays, we calculate the charmless decays 
of Ip" by virtue of the S- and ZJ-wave charmonia mix- 
ing scheme which was proposed to explain the large 



partial width and the pn puzzle. We find 



that this leads to a possible large branching fraction, up 
to 13%, of the charmless final state in tp" decays. Al- 
though the calculation is semi-quantitative, it demon- 
strates that a large charmless branching fraction in ip" 
decays can well be explained. 
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